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Abstract 

We describe a general derivation scheme for the Vlasov-type equations 
for Markov evolutions of particle systems in continuum. This scheme is 
based on a proper scaling of corresponding Markov generators and has an 
algorithmic realization in terms of related hierarchical chains of correla- 
tion functions equations. Several examples of realization of the proposed 
approach in particular models are presented. 
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1 Introduction 

Dynamical processes in many-body systems are often approximately described 
by kinetic equations, see, e.g., the excellent reviews by H.Spohn [35], [3H| . A fa- 
mous example of such equations is the Vlasov equation for a plasma, see e.g. 
[33] . |34j . The Vlasov equation in physics describes the Hamiltonian motion of 
an infinite particle system in the mean- field scaling limit, thereby taking into 
account the influence of weak long-range forces. The convergence in the Vlasov 
scaling limit was shown by W.Braun and K.Hepp [J] (for the Hamiltonian dy- 
namics) and by R.L.Dobrushin \§\ (for more general deterministic dynamical 
systems). Note that the resulting Vlasov-type equations for particle densities 
are considered in classes of finite measures (in the weak form) or integrable func- 
tions (in the strong form). The latter means, in fact, that we are restricted to 
the case of finite- volume systems or systems with zero mean density in an infinite 
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volume. A detailed analysis of Vlasov-type equations for intcgrable functions 
presented in the recent paper by V.V.Kozlov [30] . 

The main aim of this paper is to study Vlasov-type scaling for some classes 
of stochastic evolutions in continuum. Here we have in mind, first of all, spa- 
tial birth-and-death Markov processes (e.g., continuous Glauber dynamics) and 
hopping particles Markov evolutions (e.g., Kawasaki dynamics in continuum). 
Note that the approaches to the Vlasov scaling mentioned above seems to be 
quite difficult to apply to stochastic dynamics considered here (even in a finite 
volume) due to some essential reasons. For these processes, the possibility of 
their descriptions in terms of proper stochastic evolutional equations for particle 
motion is generally speaking absent. This, together with a possible variation of 
the particle number during the evolution, is an essential trouble in the applica- 
tion of the general Dobrushin's method. 

Therefore, we shall look for an alternative approach to derive the kinetic 
Vlasov-type equations from stochastic dynamics. Contrary to the classical 
derivation of the Vlasov-type kinetic equations from the Hamiltonian dynamics, 
we do not prove the law of large numbers for the corresponding processes. We do 
not even need to show the existence of the corresponding microscopic rescaled 
processes. Our main idea is to study the evolution of states (distributions) of the 
system in terms of the corresponding chain of hierarchical equations. As pointed 
out by H.Spohn [35], the correct Vlasov limit can be easily guessed from the 
BBGKY hierarchy for the Hamiltonian system. Such heuristic derivation does 
not assume the integr ability condition for the density, but until now, it could 
not be made rigorous due to the lack of detailed information about the proper- 
ties of solutions to the BBGKY hierarchy. We would like to stress that different 
classes of initial data are not only mathematical tools for the rigorous study of 
the problem. They describe different physical situations in related microscopic 
models. The zero average density systems were considered in [1] by means of 
heuristic limit transition in the corresponding hierarchical equations for corre- 
lation functions. The framework we are working in is nonzero average density 
which is related to the case of bounded correlation functions. Our approach 
is based on Spohn's observation applied in a new dynamical framework. More 
precisely, we already know that many stochastic evolutions in continuum admit 
effective descriptions in terms of hierarchical equations for correlation functions 
which generalize the BBGKY hierarchy from Hamiltonian to Markov setting, 
see, e.g., [17] and the references therein. Moreover, these hierarchical equations 
are often the only available technical tools for the construction of corresponding 
dynamics in several models [2T], [23], [II]- 

In Section 3 we propose a general scheme for the Vlasov scaling of stochastic 
dynamics for interacting particle systems in continuum. This scaling is actually 
of mean-field type which is adopted to preserve the spatial structure. Addition- 
ally, we scale the class of initial distributions at the level of the corresponding 
correlation functions. The scheme we use has also a clear interpretation in terms 
of scaled Markov generators. An application of the considered scaling leads to 
the limiting hierarchy which possesses a chaos preservation property. Namely, if 
we start from a Poissonian (non-homogeneous) initial state of the system, then 
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this property will be preserved during the time evolution. The main observa- 
tion which appears at this point is the following. A special structure of the 
interaction in the resulting virtual Vlasov system gives a non-linear evolutional 
equation for the density of the evolving Poisson state. It is for the first time that 
macroscopic Vlasov-type equations are obtained from the microscopic infinite- 
particle systems in an unbounded region of non-zero average density using the 
corresponding system of hierarchical equations. 

Section 4 is devoted to the application of the general scheme to a wide 
class of birth-and-death and hopping particles processes. We state conditions 
on structural coefficients in the corresponding Markov generators which give 
a weak convergence of the rescaled generators to the limiting generators of the 
related Vlasov hierarchies. As a result, we may compute the limiting Vlasov- 
type equations for the considered processes leaving the question about the strong 
convergence of the hierarchy solutions open. In Section 5 we present a collec- 
tion of particular examples of the resulting Vlasov equations for several concrete 
models. Note that each of the examples considered creates its own non-linear 
equation for the density in the discussed scaling. These equations include con- 
volution operators as a common point of their structure. To our knowledge, any 
general results concerning properties of solutions to such kind of non-linear evo- 
lutional equation are absent. This is an exiting mathematical problem strongly 
motivated by concrete models of interacting particle dynamics. 

Many problems of the (mathematical) population biology concerns interac- 
tions between populations of different types. Our technic, in fact, covers this 
case. In particular, one can derive spatially inhomogeneous non-linear equations 
of the Lotka-Volterra type in the Vlasov-type scaling. On the other hand, one 
may apply this approach to the so-called continuous Ising model (Potts model) . 
We explain these results in forthcoming papers [TU] , [13] • 

Note that control of convergence of the Vlasov scalings for the solutions to 
considered hierarchies is a difficult technical problem which shall be analyzed 
for every particular model separately. Our results in this direction concern 
two classes of models: Glauber dynamics in continuum and a spatial ecological 
model (so-called Bolker-Dieckmann-Law-Pacala model). Due to their techni- 
cally complicated character, these results will be published in separated works 

2 Basic facts and notation 

Let B(R d ) be the family of all Borel sets in R d , d > 1; B h (R d ) denotes the 
system of all bounded sets in B(R d ). 

The configuration space over space M. d consists of all locally finite subsets 
(configurations) of M. d , namely, 

T = r Rrf := { 7 C R d | | 7 n A| < oo, for all A e B h {R d )}. (2.1) 

The space T is equipped with the vague topology, i.e., the minimal topology for 
which all mappings r 3 7 1— > X^e-/ f( x ) £ K are continuous for any continuous 
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function / on R d with compact support; note that the summation in Ylxe-y f(x) 
is taken over finitely many points of 7 which belong to the support of /. In 
[2"U] . it was shown that T with the vague topology may be metrizable and be- 
comes a Polish space (i.e., complete separable metric space). Corresponding to 
this topology, the Borel cr-algebra B(T) is the smallest er-algebra for which all 
mappings T 3 7 h-> G N := N U {0} are measurable for any A G Sb(M d ). 
Here 7a := 7 PI A, and | • | means the cardinality of a finite set. 

The space of n-point configurations in an arbitrary Y G B(R d ) is defined by 

r^ l) := {n C Y I \ V \ =n}, n G N. 

We set also TyP :— {0}. As a set, Ty^ may be identified with the symmetrization 
of _ 

y = {(an, .. .,x n ) G Y n \ x k ^ Xl ifk^l}. 

Hence, one can introduce the corresponding Borel cr-algebra, which we de- 
note by S(r^). The space of finite configurations in an arbitrary Y G B{W l ) 
is defined by 

T . Y := □ r< B) . 

nGN 

This space is equipped with the topology of disjoint unions. Therefore, one can 
introduce the corresponding Borel cr-algebra B(To t Y)- In the case of Y = M. d we 
will omit the index Y in the notation, namely, Tq ■= F ^d, :— rfc . 

The restriction of the Lebesgue product measure (da;)™ to (T^ n \ fi(r( n ))) we 
denote by rrv- 71 '. We set :— S^y The Lebesgue-Poisson measure A on r 
is defined by 

00 

A: =E- m(n) - ( 2 - 2 ) 

For any A G Bh(R d ) the restriction of A to Ta := ^q,a will be also denoted 
by A. The space (r,£>(r)) is the projective limit of the family of spaces 
{(Pa, S(rA))} Ag g b ^ R(i ^. The Poisson measure tt on (r,6(r)) is given as the 
projective limit of the family of measures {7r A }AeB b (R d ): where ir A := e~ m ^X 
is the probability measure on (Fa,£>(Fa)). Here m(A) is the Lebesgue measure 
of A g 23 b (R d ). 

For any measurable function / : R d — > M we define a Lebesgue-Poisson 
exponent 

ex (f,v)-= II ^ er °; e A (/,0):=l. (2.3) 

x£Lr] 

Then, by (E2|), for / G L 1 (K d , dx) we obtain e\(f) G L^ro.rfA) and 

e x (f,r))d\(ri) = exp{ / f(x)dx). (2.4) 
r UR d J 
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A set M £ B(T ) is called bounded if there exists A £ B h (R d ) and N £ N 
such that M C U^Lo^A • The set of bounded measurable functions with 
bounded support we denote by Bb s (T ), i.e., G £ Bbs(ro) if G \r \M = f° r 
some bounded M £ B(T ). Any £>(r )-measurable function G on To, in fact, is a 
sequence of functions } nE jj where is a S(r'™')-measurable function 

on r'"l We consider also the set J-" cy i(r) of cylinder functions on I\ Each 
.F £ Tcy^T) is characterized by the following relation: F(j) = F \r A (7a) for 
some A £ B h (R d ). 

There is the following mapping from _Bbs(ro) into J-" cy i(r), which plays the 
key role in our further considerations: 

KG( 7 ):=J2G(V), 7eT, (2.5) 

where G £ Bb s (T ), see, e.g., [HJ [3T] [35] ■ The summation in (I2.5[) is taken 
over all finite subconfigurations r\ £ To of the (infinite) configuration 7 £ T; 
we denote this by the symbol, r\ <g 7. The mapping X is linear, positivity 
preserving, and invertible, with 

K- x F{j ] )x=YJr-^ m ni\ ver a . (2.6) 

Here and in the sequel inclusions like £ C 77 hold for £ = as well as for £ = rj. 
We denote the restriction of K onto functions on To by i^o- 

For any fixed C > 1 we consider the following Banach space of B(Tq)- 
measurable functions 

' \\G\\ C := [ \G(r,)\CMdX( V ) < 00). (2.7) 



G:IW 

A measure /i £ Mj m (T) is called locally absolutely continuous with respect 
to (w.r.t. for short) the Poisson measure tt if for any A £ Bb(R d ) the projection 
of fi onto Ta is absolutely continuous w.r.t. the projection of tt onto T\. By 
[19], in this case, there exists a correlation functional k^ : To — > K+ such that 
for any G £ i?bs(ro) the following equality holds 



(A-G)( 7 )d/*(7) = / G{rf)kMdm- (2-8) 
r Jr 

The restrictions of this functional on Iq , n 6 No are called correlation 
functions of the measure [i. Note that fc^ "* = 1. 

We recall now without a proof the special case of the well-known technical 
lemma (cf., |28] ) which plays very important role in our calculations. 

Lemma 2.1. For any measurable function H : To x To x To — > R 

J '£H(Z,r 1 \Z,r 1 )d\(r 1 )= f f H (£, V , r? U dX (£) dX (77) (2.9) 

onZy if both sides of the equality make sense. 
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3 General scheme 



In this section we introduce the notion of the Vlasov scaling for Markov dynam- 
ics of IPS on configuration spaces. 

We assume that our system evolves in time due to some mechanism whose 
details will be specified for concrete models. Suppose that the initial distribution 
of particles in our system is a measure fig £ .M.j m (T), with correlation function 
fco- Let fit S A4 1 (r) be the distribution of particles at time t > and kt be its 
correlation function. One should note that if evolution (fJ-t)t>o is ruled by an a 
priori given Markov process on T (i.e. if such a Markov process exists), then fj, t 
is a solution to the following Kolmogorov equation: 



where L* is the operator adjoint to the generator of functional evolution, i.e., 



Of course, one should be careful about the functional and measure spaces to 
have all the above-introduced operators properly defined. We postpone careful 
definitions of all these objects until the introduction of concrete models. 

Now, assume that the evolution of correlation functions (k t )t>o, correspond- 
ing to (nt)t>o, first of all exists, and is the solution of the following evolutional 
equation 



where L is the generator of a semigroup T t on some functional space which 
includes all bounded functions (or bounded with some weight) almost every- 
where (a.e.) w.r.t. the Lebesgue-Poisson measure A. In many applications this 
space may be taken to be Kc '•= {k : To — >• K | k ■ C^' 17 ' <E L°°(A)} for some 
fixed C > 1. Let us stress that (|3.1j) is nothing else but a hierarchical system 
of equations corresponding to the Markov generator considered. This system 
has the same meaning as the BBGKY hierarchy in the case of Hamiltonian 
dynamics. 

The first important step on the way to construct the Vlasov scaling concerns 
the proper rescaling of the initial state of the system. Or, equivalently, in the 
language of correlation functions it means the proper rescaling of the initial 
conditions of (|3.1I) . 

More precisely, at the beginning we rescale fco with parameter e > in such 
a way that the resulting functions kg as e — > behave as follows: 






(3.1) 



^0,rcn( 7 ?) : ~ £^k { 



('?) -> r (v), £ 



-> 0, V e r , 



(3.2) 



G 



where the function tq is a subject of choice for concrete examples and aims. In 
general, it has to be a bounded function also (or bounded with some weight) 
a.e. w.r.t. the Lebesgue-Poisson measure. 

Remark 3.1. In the case of r (r]) = e\(po,r]), r\ s r , po ■ M d — > (0,+oo) the 
assumption about the rescaling of the initial condition means heuristically the 
following: Pg\ cn —> 7r po , where p^\ cn has a correlation function e^l/c^^). 

It is clear that such a rescaling of the initial solution for (|3.ip leads to a 
singular function w.r.t. e > 0. In applications, this fact can be interpreted as 
the growth of density of the system with e — > 0. 

We have to consider (and it is our second step) some proper scaling of the 
generator in (|3 .If) : 

L A i— ► L A . (3.3) 

The concrete type of this scaling will depend on L A . In the next sections we 
consider several types of generators and corresponding scalings. Suppose that 
there exists a solution of the functional evolution 



dk. 



*— - LH[ e) 



dt 6 1 (3.4) 

tOOl _i.00 
> K t \ t =o ~ K o 

We expect (and this will be shown in the concrete models for the concrete 
scalings in forthcoming papers) that this solution will be also singular w.r.t. 
e > 0, hence, this solutions will be in functional spaces depending on e. 

Moreover, we should choose the type of scaling (|3.3|) which guarantees that 



the order of this singularity will be the same for the initial function k^ . Namely 
(and it is our third step on the way to realize the Vlasov scaling) we consider, 
cf. (53) . 

^L(v)--=£ l M E) (v), ier„, (3.5) 

and we want to show that 

fcfLW ^o,,,er„. (3.6) 

In fact, (|3.5[) means that we consider a renormalized version of the evolution 
equation (|3.4p : 

aK _t 1 icn _ T A ,(e) 

dt e ' ron 4 ' ren (3.7) 



JL(=) I _!.(=) 
■ A 't,ren| t= o — O.ren 



where 



L*= E ML?e-W (3.8) 



Therefore, informally, we want to show that the solution of the evolution 
equation (|3.7j) converges (in a proper sense) to some function rt which satisfies 
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the Vlasov hierarchy 

< dt (3.9) 

Recall again that the choice of the scaling (|3.3|) is prescribed by the model. 
Having applications in mind, it is important to consider the case of ro(rj) = 
Ca(pO) 7 ?) and the scaling (|3.3I) which leads to r t of the same type, i.e., 

n(rf) = e\(Pt,r)), V^^o- 

The latter means the so-called chaos preservation property of the Vlasov hier- 
archy. Equation (13.91) in this case implies, in general, a non-linear equation for 
Pt- 

^Pt(x) = v(jH)(x), xeR d , (3.10) 

which we will call the Vlasov-type equation. 

To describe this scheme in a more analytical way, we use the language of 
semigroups. Suppose that we know the mechanism of the evolution of our system 
given by the Markov pre-generator L. Let L be defined at least on functions from 
•^cyi (r) and L = K^ 1 LK be the corresponding descend mapping on functions 
from £?bs (ro)- Let us fix the duality between functions on Tq 

({G,k)) = [ G(rj)k( V )d\( V ), (3.11) 

and consider the mapping L being the dual to L w.r.t. (|3.11l) . 

Assume that L can be extended to a generator L. We want to construct a 
scaling of the generator L, say, L £ , e > 0, such that the scheme described above 
will be covered. Assume that we have a semigroup T £ (t) with a generator L £ = 
K^ 1 L £ K in some functional space over Tq. Consider the dual semigroup T A (t) 
which corresponds (in a proper sense) to L A . As we said before, we consider an 
initial condition of (|3.4|) with a singularity in e, namely, (ff) ~ e _ ' r? 'ro(?j) ) 
e — > 0, 77 G Tq with some function ro, independent of e. First of all, we have 
to choose such a scaling L n> L e for which T^(t) preserves the order of the 
singularity: 

(T £ A (t)k a E) )( V )^£^r t ( V ), £^0, ryeLo. (3.12) 

And the most important is that the dynamics ro H> r t should preserve the 
Lebesgue-Poisson exponents: if ^0(77) = ex(po,rf) then rt(rj) = e\(pt,r/), where 
pt is satisfied (|3.10|) . 

Now let us close our construction with the evolution of states in this scheme. 
Let us consider for any e > the following mapping of functions on Tq 

(R e r)(ri):=£^r(v)- (3-13) 

This mapping is "self-dual" w.r.t. duality p. lip , moreover, i?" 1 = R E -i. Then 
we have k { Q s) - R B -ir , and we need r t - R £ T e A (t)k { Q e) - R £ T £ A (t)R £ -ir Q . 
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Therefore, we have to show that for any t > the operator family R £ T A (t)R £ -i , 
e > has a limiting (in a proper sense) operator U(t) as e — > and 

U(t)e x (p ) = e x { Pt ). (3.14) 

But, informally, T A (t) = exp{tL A } and R £ T A (t)R £ -i = exp{tR £ L A R £ -i}. In 
fact, we need the existence of an operator V A such that exp {tR £ L A R £ -i } — > 
exp {tV A } =: U(t) for which (|3.14p holds. Therefore, a heuristic way to produce 
the scaling L t— y L £ is to demand that 

l^l {§t e ^ Pt,v ^ ~ L t^n e x(Pt,v)J = 0, V e T (3.15) 

if p t is satisfied (I3.10[) . The point- wise limit of L £rcn will be the natural candi- 
date for V A . 

Sometimes, to show convergence of solutions of evolutional equations in some 
functional spaces it is much simpler to work with the operators L £ , lcn and V 
which are pre-dual to L A ren and V A w.r.t. the duality (|3 . 11 [) . Note that (|3.8j) 
implies 

^e.ren = R £ -^L £ R £ , (3.16) 

and V should be the point-wise limit of -L E)ren - 

4 Generators of birth, death, and hopping 

Through out this section we consider generators of two types for continuous 
models: the birth-and-death generator Lbad = L~ + L + and the hopping gen- 
erator Lhop, where for any F 6 J-" cy i(r) 



(L-F) ( 7 ) := ]T d (x, 7 \ x) [F( 7 \x)-F (7)] , (4.1) 

(L+F) (7) := f b (x, 7) [F (7 U x) - F (7)] da, (4.2) 

(L hop F) (7) := [ c (x, y, 1 )[F( 1 \xUy)-F (7)] dy. (4.3) 
xe-y jRd 

Here b, d, c are measurable functions of their variables and, additionally, b and c 
are locally integrable function of the first and second variables, correspondingly. 
These conditions guarantee that (|4.1l) - (|4.3p are well-defined on .F cy i(r) since for 
any F 6 J" cy i(r) there exists some A e B h {R d ) such that F(j \ x) = F^) for 
any x £ jA", F(j Ui) = ^(7) for any x £ A c , and F(j \ xUy) = Ffa) for any 
x G 7A=, y e A c ; as result the sums in (|4.1j) and (|4.3j) are over finite set 7a and 
the integrals in (|4.2j) and (|4.3p are over bounded set A. 

We may denote L~ — L~(d), L + = L + (b), Lh op = Lh op (c). Assume that we 
have some scaling of rates b,d,c, say, b £ ,d £ ,c e , correspondingly; e > 0. Then, 
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let us consider the following scaling of Lbad and Lh op : 

L bad , e = L-(d e ) + e- 1 L+(b e ), (4.4) 
Lho P ,e = L hop (c e ). (4.5) 

Remark 4.1. In a conservative system with a generator like (|4.3I) which preserves 
the "number of particles" during an evolution the Vlasov-type scaling usually 
means decreasing of the intensity of the interactions between elements of a 
system together with increasing of correlations in the initial state. However, in 
a non-conservative birth-and-death dynamics with a generator L D ad we need an 
additional increasing of the birth intensity to preserve the influence of the birth 
part in the limiting Vlasov hierarchy. Note that the necessity of the concrete 
factor e _1 in (|4.4p is clear a posteriori only (see Proposition I4.5|) . 

Suppose that there exists three families of measurable functions on IV D x e , 
B { x e \ Ci £ ) y , £ > 0, {x, y} C R d , such that 

4 (a; l7 ) = (KD^)( 7 ), b E (x, 7 ) = (i^i £) )(7), 

c £ (x,y,j) = (KC£i)( 1 \x). 

Note that, in general, C XtV ^ C VtX . 

Proposition 4.2. The following formulas hold for any k € i?bs(ro) 
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Then, gU) and gTJ follow directly from flS^). □ 
Let pt, t > be measurable functions on M. d . The explicit formula 

^e A (p t , rj) = e A (Pt, V \ x) -T^pt 0) (4.8) 

together with our "demand" (13. 15)) induce us to state the following corollary. 
Corollary 4.3. Let p be a measurable function on M. d . Then 

(Lt^, m MP)) (V) (4-9) 
= -J2e x (p,r 1 \x)p(x) f e x (p,0 E ^ U dX (0 

+ J2e x (p,r 1 \x) f e x (p,0 £ er-l^'^? (a, U f) dA (0 ; 



(^Vp.^xO,)) (v) (4.10) 
= J2e x (p,r)\x) [ p(y) [ e A (p,0 V £ H?I C$ (w U dA (0 

-J2e x (p,r 1 \x)p(x) [ [ e x (p,0 £ e-KIC^U^dA 

Proposition 4.4. Suppose that for any {x,y} C R d , C To 

3 lim E e-^D^ (« U £) = ltae"^!^ (0 =:^(0. ( 4 - n ) 

3 lim 2 ^'4 E) U = lim e-l?!^) (£) =:i#(£), (4.12) 



3 ]™o 51 *~®<%l ( w u o = H £ ' m °^ (e) = :C ^v(0- ( 4 - 13 ) 



Then, our "demand" (13. 1 5[) holds. More precisely, 

(V h A ad k) (v) ■= Jim (L^^k) (V) (4-14) 



2*(CU(Tj\a:))Sj(e)^(0. 

x£rj 
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and if p t is the solution of the equation (|3.10p with 

v{ P )(x) = v had (p)(x) =-p{x) [ e x (p, D v x (£) d\ (0 

e A (p,0^(0dA(0. ( 4 - 15 ) 



Q 

then, -^-e\{p t ,n) = (V haA e\(p t )) (rj). Analogously, 



01 

(V h A op k) („) := lim (Lt p , e , rcn fc) (V) (4.16) 
= E / / k(tU(r ) \x)Uy)CZ x (0d\(0dy 

- W fc(eur?) / C^(0d»dA(0, 
xen Jr Jn* 

and if p t is the solution of the equation (|3.10p with 

v(p)(x) = v hop (p)(x)= [ p(y) [ e x (p,0CY x (0 d\(0dy (4.17) 



-p(x) e A (p,0 / 0^(0^(0, 
•/r ^R d 

<9 

#*en, T^e x {p t ,v) = (V hop e x (p t )) (rj). 

Proof. The equalities (I4.14p and (|4.16l) are direct consequences of the Proposi- 
tion gTH and the conditions ([4TTT]I - ([4TT5]) . Taking the limit in (|4"^| and (|4~TU)) 
as e — >• and using (I4.8[) we obtain the statement. □ 

And now we present the explicit expressions for the corresponding operators 
L Ej ren and V. 

Proposition 4.5. For any G 6 B^Tq) the following formulas hold 

(£ ba( Wc„G)( ?7 )=-]r J2 G(ZU X )J2z- l(vVMl Di S H"U(r ] \x)\0 



+ V / G(^ux) V E -to\*l£W( w u»AO<te; (4-18) 
(W, rc „G) (??)=£ Yl I G^Ux)^2e-^y^C^ x (u;U(r 1 \y)\0dx 
"E G K)E E J j-^CMiuUTiWdv. (4.19) 



£C?7 x6£wC£\a: 



12 



//, additionally, (|4.11l) - (|4.13p hold, then, 

(FbadG) fo) = 53 G (0 E ^ fa \ o 



+ 53 / G^U^B^U)^; (4.20) 
(VW?)fa)=E E / G^U^G^^U)^ 

"E G ^E / CZ y { V \Z)dy. (4.21) 



£07 xef 

Proof. We may obtain these formulas directly from the duality p. lip and the Lemma l2~Tl 
Namely, for any G £ £?bs(ro) we have 

G (T)) (ibad, e,rcnfc) (v) ^A (v) 

To 

-f G(V) [ fc(£U77)E E e- |fl £>i e) (wU0dA(0dA(»j) 
+ [ G( V ) I 52k(tU(ri\x)) e-MBV(uUt)d\(e)d\(ri) 



mG(r 1 Ux)k^Ur 1 Ux)y] e~^D^ (u U £) dirdA (£) dA 

wO? 

/ / / G(r 1 Ux)k(ZU V )Y,£- m Bi £) ("U0dxd\(0d\(r 1 ) 
Jr Jr Jul* wCr? 

/ / E G ^ Ua; ) fc ^ Ua; ) E 5 ''^'- ^ (w U £ \ 77) dardA (0 
•/r •/*« „ cC wCr) 

/ ^ / G (77 U *) fc (0 ^ e- 1« \"l ( W U e \ t,) dzdA (0 



/EE G(ryUx)fc(?)E eH ^ V ' l ^ ) ( wU ^\ ;c \^ dA «) 
/ 53 / G(r)U x)k (£) 53 £- |?V7l 4 e) (w U £ \ ft) dxdX (0 , 



which implies (|4~T8|) . To get (|4~20f we may proceed in the same way or just let 
£ ->• in (|4~T8|) . Then (|+Tll - (l4T2l) together with equality 

53 53 G(tUx)DV((r 1 \x)\0 = J2 G (Oj2 D x( r l\0 
provide P~2"0"|) . 
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Analogously, for any G G Bbs(ro) we have 
k) ( V )d\(r,) 



[ G(r,Ux) [ f k^\Jf 1 Uy)J2e-^C^%(u J u0dX(0dxdydX(r 1 ) 
/ G( V ) [ k(tU V )J2J2 I e-^C^(u;UOdyd\(Od\( V ) 



I I I Y, G ^ Ux ) k ^y)Y, £ ^ C vl { - uj{J t\ri)dxdyd\{0 
Jr Jr« Jn* nCi ^ 

-f E G ^) fc ^E E / e-^CViuU^dydXiO 

I k & I EE G( f) u a: )^ £ -i^i^( W ue\!/\)i)«A(e) 

-/ fc(C)E G ^E E / e-^ v{ Cii(uUt\v)dyd\(0, 



which implies (|4.19[) . To get (|4.21l) we may proceed again in the same way or 
just let e ^ in (|4~19l) and use (|4~13l) . □ 

In the next Section we consider concrete examples for the operator L. 



5 Examples 

As we have seen in the previous section, the sufficient conditions (|4.1ip - ||4.13p 
have identical structure for death, birth and hopping parts. Therefore, to 
present explicit expressions for L^ len , V A and others we may proceed in the 
following manner. Let 0(7) = [KA)(~{), where A is a measurable function on 
To; let a £ = KA e be some scaling of a and A, e > 0. Below we consider different 
types of the function a (linear, exponential etc.) and present possible scalings 
such that for any {ry, £} C T 

3 lim (wUf) = lim e~^A e (£) =: A v (£_). (5.1) 

And after that we may apply this results to the our situation when A £ depends 
additionally on x, y £ R d . 

1. Let 0(7) Eael. Then A(jf) = a ■ O'' 7 ' and we don't need scaling at all: 
if a e = a then A e {rj) = a ■ and dHTTJ) holds with A v (£) = a ■ 0^. 
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2. Let 0(7) = £ f(x) with some / : R d ^ K. Then Afa) = X {j)=M} /(z). 
We consider the scaling / M> e/ for which 



Urn e ~^ A e (w U = lim ^ £" kl X {w ue={x}}£/(^) 
= lim £" l?l X{£={z}}£/(z) + lim ^ ef(x) = X{ i= {x}}f(x) =: A v (0- 



3. Let 0(7) = exp{ £ f(x)}, f : R d ^ K. Then = e A (e / - 1,77). We 
consider the same scaling / i-> e/ for which 

lim V £-l ? U £ (w U = lim V £- |?l e A (e ef -l.wUf) 

= lime A lime A (e^,r?) = e A (/,0 =: 

4. Let 0(7) = £ gi x iV) f° r some (non-symmetric, in general) function 
5 on R d x R d . Then 

^ (»?) = XI ^o -1 (X! ^ y fa \ x ) = H X{t,\x={j/}}5 (a;, y) 

=X{r,={z,y}} [5 y) + 5 (y, a;)] . 
We consider the scaling g i-> e 2 g. Then 



lim ^£-l«U e (wUO 

lim ^ e _l?l X{ w ue={x,y}}£ 2 [<7 0, y) + 5 (2/, x )] 



lime 2 X{i={ a: , ? /}}£ 2 [g(x,y) + g (y, x) 



e->0 



=X{t={x,y}} [5 (a:, 2/) + 9 {v, x)] =: A v (£). 
5. Let 0(7) - E /(^)cxp{ £ fffotf)}, / :M^l,.?:l d xE d 4 

2:67 yE~/\x 

Then 

A(r ? ) = ^/( a; )e A (e^-)-l,r ? \ a; ). 
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Let us consider the scaling / i-> ef, g n> eg. Then 
lim^e-l«U £ (wUO 

= lim^£-l«l ^ £ /(j;)e A (e*'»-l, ll ;U^i) 

e £g(x,-) _ j 



iToEE^ ( e£9(x '' } -l,u\x)ex f - ,£ 

+ ]™ E 51 (e £9(x '° - 1, w) e A ( 



e eg(x,-) _ ^ 



e sg(x,-) _ i 



= iTo E (e efl(x, °. »7) eA I 

6. Let 0(7) = ( £ /(*)) CX P { E «?(»)}, /,S = R d ^ R- Then 

x€7 2/67 

A (77) = (x { .= {x} }/(a;) * e x (e 9 - 1, •)) (v) 

= ]T /(x)e A (e» - 1, r,) + E (e 9 - 1, r; \ x) . 

Let us consider the scaling / i-> £,f, 5 n> £<?. Then 
lim^£-l«U £ (wUO 



= lim^£- | « l £ #h(e £s -l,wU() 
+ li m y£-l«l V £/(x)eA(e £S -l,a;Ue\a;) 

= E E s f( x ) e * - L w ) ^ > 

— ^— >o 

+]™ EE £ /( x ) e ^ e£9 - 1 ' cj \ a; ) e * ( £ -^^ 

— 

. . , , , - 



: E/(^ (gA\x)=--A v (0. 
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Now we consider different types of birth-and-dcath and hopping models with 
rates which have one of the forms considered above. Using explicit expressions 
for A and scaling for each concrete model we have the expression for A e (which 
is or Cx)y) and may easily obtain expressions for £ e . r0 n and L^ mn 

from (|4.6[) or (|4.7j) . Using expressions for A v (which is D^, B^ or C^ y ) we 
may obtain expression for V and U A as well as the form of v also from the 
Propositions 14.41 and 14.51 Let us turn to these concrete examples. We present 
the Vlasov-type equations only. 

Example 5.1 (Surgailis model). This birth-and-death model describes indepen- 
dent appearing and disappearing points from a configuration after exponentially 
distributed random times. The corresponding dynamics was considered in |37] . 
[55] ; the generator may be given for F G J" cy i(r) 

(LF)( 7 )= m y2[F( 1 \x)-F( 1 )} + a f [F( 7 U x) - F(j)]dx. 

The scaling m H ► m, a h- > e~ x o leads us to the following Vlasov-type linear 
equation 

® t \ t \ 

-Q t Pt{x) = -mpt(x) + a. 

Example 5.2 (Contact model). This model was considered in [5^ (for further 
investigations see [22], [II])- The model describes independent death of the 
members of a configuration, and, on the other hand, production of new mem- 
bers of the configuration by the existing ones. This is the simplest model for 
ecological population dynamics. Note that a similar model was considered al- 
ready in [7] as a particular case of a spatial branching process in continuum. 
The generator is given on T cy i(T) by the expression 

(LF)(7) =m'£[F( 1 \x)-F( 7 )} 

+ A V f a(x-y)[F( 7 Uy)- F( 7 )]dy. 

The described scaling m i— > m, X h-> e _1 A, a M> ea (that means that L = L e ) 
provides the linear Vlasov-type equation also 

d 

— p t (x) = -mpt(x) + X(p t * a)(x). 
Here and below * denotes the usual convolution in M. d . 

Example 5.3 (Social model). This model was considered in [5]. It describes 
birth-and-death process with migration from some "reservoir" and competition 
between members of a configuration. The generator is given on J- cy \(T) by the 
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expression 

(^)(7) = £ E o(»-v)[J , (7\«)--F'(7)] 

2:67 y£~t\x 

+ a [F{jUx) - F(j)]dx. 

The described scaling a 1— > ea, er i— > £ <r provides the non- linear Vlasov-type 
equation: 

d 

glPt(x) = -pt{x){p t * a)(x) + a. 

Exam/pie 5.4 (Bolker-Dieckmann-Law-Pacala model). This model of popula- 
tion ecology was considered in [2J, [3], [5]. Rigorous mathematical studying of 
this model was done in |14j . The individual of a population may die indepen- 
dently as well as due to competition for resources; any individual may produce 
a new one also. The generator is given on J-" cy i(r) by the expression 

(lf)( 7 ) = ]T( to + E a 'i x -y))l F (y\x)-F('y)} 

x£y ye-/\x 

+ A E/ a + (x~y)[F( 1 Uy)-F( 1 )]dy. 

The scaling >-> ea il 1 m 1— > m, A 1— > e~ 1 A gives the following non-linear 
Vlasov-type equation: 

d 

■^lPt( x ) = - m Pt( x ) ~ Pt{x){Pt * a )(x) + X(p t * a + )(x). 
at 

Note that in the space-homogeneous case we obtain the logistic-type equation 
^p t = (A(a+) - m - {a~}p t )pt, 

where (a ) = J R d a [x)dx. For a rigorous proof of convergence in this scaling 
see [TT]. 

Example 5.5 (Contact model with establishment). In this model the above de- 
scribed contact dynamics is improved by taking into account the depressive role 
of the establishment. Namely, the probability for a newborn member to survive 
in a new place is smaller if there are more particles near this new place. In the 
language of a generator we describe this by the following expression 

(LF)( 1 )=mJ2[F(j\x)-F(j)} 

ai£7 

+ A E / y)e~ [F( 7 U y) - F(~f))dy. 
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The scaling m h- > to, A i— > e A, a n- ea, i— » provides the following non- 
linear Vlasov-type equation 

^pt(x) = -m Pt (x) + A(a * p t )(^)e- (0 * pt)(x) . 

Example 5.6 (Contact model with fecundity). This model describes influence of 
competition for resources on birth intensity. Namely, if there are many existing 
members near a "parent" , the probability to sent offspring for it is smaller. We 
consider the following expression for the generator 

(LF)( 7 )=m^[F( 7 \x)-F( 7 )] 

+ A^V^e-x*^^ f a(x-y)[F( 7 Uy)-F('y)]dy. 

The previous scaling to h- > to, A i— > e _1 A, a t— > ea, i— > e^> yields another 
non-linear Vlasov-type equations 

JU(x) = -m ft (x) + A(a * (p t e-^)))(x). 

OT 

Example 5.7 (Dieckmann-Law model). This model, as well as the model from 
Example 15.41 describes ecological population evolution. However, appearing of 
new offsprings is proportional to the number of existing members of a popula- 
tion. The generator is given by the following expression 

(£F)(7) = £(m+ E a-(x-y))[F(j\x)-FW] 

x<£~/ y£j\x 

; f a+(x-y)(x + V b(x-u))[F( 7 Uy)-F( 1 )]dy. 



+E 

Note that without competition (a - = 0) this model explodes, namely, the mean 
value of the number of members in any bounded region becomes infinite after 
finite time; otherwise, if the competition kernel a~ is "stronger" than the kernel 
6 this effect is absent (for details see [5]). After scaling a^ 1 H> ea^, to h-> to, 
b H» eb and 1 H» (before the whole birth term) we obtain the following 
non-linear Vlasov-type equation 

d 

TnPt{x) = -mp t (x) - pt(x)(p t * a )(x) + X(p t * a + )(x) + (((& * p t )pt) * a + )(x). 
at 

Example 5.8 (Glauber G + dynamics). This model is a continuous analog of the 
Glauber dynamics on a lattice. It was considered in a couple of works, see, e.g., 
[21], E3], EH], US, [IS]- The generator of this model is given by 

(LF)( 7 ) = y2[F( 7 \x)-F( 7 )} 



x£7 



e-y [F( 7 U j/) -F(j)]dy. 
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Here z > is an activity parameter and is a pair potential. This generator 
has a reversible measure, namely, the Gibbs measure with parameters z and cj> 
(see, e.g., [23], [TB] for details). The scaling m h-> m, z M> £ _1 z, <f) ^ e<f) yields 
the following non-linear Vlasov-type equation 

^- t p t (x) = - Pt (x) + ze-^*^\ 

For a rigorous proof of the convergence in this scaling see [12] . 

Example 5.9 (Glauber G~ dynamics). This model is similar to the previous one, 
see, e.g., [UJ, [26]. 

(^)(7) - E eE " £7 0(a; - u) ^(7 \ x) - F( 7 )] 
+ z / [F( 7 U y) - F( 7 )]dj/. 



The same scaling as before yields the similar non-linear Vlasov-type equation 
^-p t (x) = -p t (x)e^*^ + z. 

Example 5.10 (Free Kawasaki). This simplest exactly solvable hopping model 
was considered in [27]. It describes independent jumps of particles in the system. 
The generator is the following 

(LF) ( 7 ) = V / a (x - y) [F ( 7 \ x U y) - F ( 7 )] dy. 

We do not need scaling at all to obtain the linear Vlasov-type equation 
d 

—p t (x) = (p t * a) (x) - (a) p t (x) . 

Example 5.11 (Density dependent Kawasaki). In this model the intensity of a 
jump is linearly proportional to the existing population. The generator is given 
on .F C yi(r) by the expression 

(LF)( 7 )=V f a(x-y)J2b(x,y,u)[F( 1 \xUy)-F( 1 )}dy. 



The scaling b h-> eb provides the following non-linear Vlasov-type equation 
d 

-p t (x)l a(x-y) I p t (u)b(x,y,u)dudy. 



^lPt( x )^( Pt{y)a{x~y) I pt{u)b(y,x,u)dudy 
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In particular, if b(x, y, u) = b(x — u) then 



■rr.Pt {%) = {(pt{pt * b)) * a)(x) - (a)p t (x) (p t * b)(x). 



If b(x, y, u) = b(y — u) then 



— p t {x) = (p t * b) (x) (p t * a) (x) - p t (x) (p t *a* b){x). 



Example 5.12 (Gibbs-Kawasaki) . This hopping particles model was considered, 
e.g., in |26j . The generator is given by the expression 



It has a family of reversible Gibbs measures with the potential (f> and any activity 
z > 0. The scaling (j> h-> e<f> gives the non-linear Vlasov-type equation of the 
form 

d 

—pt (x) = (p t * a) (x) exp {- (p t * ((>) {x)} - p t [x) (a * exp {-p * (j)}) (x) . 

Example 5.13. In the last example we consider another type of dynamics. Let 
L describe the generator of the diffusion dynamics (see, e.g., [25], [H]), namely, 
for any smooth cylindrical function 



where A is a classical Laplace operator in K d and V is a gradient in R d . Our 
approach covers this case also. It can be shown that the scaling (j> \-t e<\> provides 
the following non-linear partial differential Vlasov-type equation 
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